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Abstract 

In this paper we develop constructive invertibility conditions for the twisted convolution. Our approach 
is based on splitting the twisted convolution with rational parameters into a finite number of weighted 
convolutions, which can be interpreted as another twisted convolution on a finite cyclic group. In analogy 
with the twisted convolution of finite discrete signals, we derive an anti-homomorphism between the 
sequence space and a suitable matrix algebra which preserves the algebraic structure. In this way, the 
problem reduces to the analysis of finite matrices whose entries are sequences supported on corresponding 
cosets. The invertibility condition then follows from Cramer's rule and Wiener's lemma for this special 
class of matrices. The problem results from a well known approach of studying the invertibility properties 
of the Gabor frame operator in the rational case. The presented approach gives further insights into 
Gabor frames. In particular, it can be applied for both the continuous (on M"^) and the finite discrete 
setting. In the latter case, we obtain algorithmic schemes for directly computing the inverse of Gabor 
frame-type matrices equivalent to those known in the literature. 



1 Introduction 

Twisted convolution arises naturally in the context of time frequency operators, more specifically in the 

treatment of Gabor frames ^|S]. The study of inversion schemes of twisted convolution has, therefore, a 

*Dept. of Electrical Engineering, Technion-Israel Institute of Technology, 32000 Haifa, Israel. 

Tel.: +972-4-8293256, Fax.: 4-972-4-8295757, Email: yonina@ee.technion.ac.il 

^Faculty of Math., University of Vienna, 1090 Vienna, Austria. Tel.: -1-43-1-4277-50693, Fax.: -1-43-1-4277-50690, 

Email: ewa.matusiakOunivie . ac . at, tobias . wertherOunivie . ac . at. 



1 



major impact on the analysis of Gabor frames. Our method is originated by the Janssen representation 
of Gabor frame operators 11 and simphfies the approach given in ^^l- A different, however, equivalent 
method for studying Gabor frame operators is the well known Zibulski-Zeevi representation ^1] based on a 
generalized Zak-transform. 

In contrast to the standard convolution, the twisted convolution is not commutative. This is opposed 
to the possibility of applying powerful tools from harmonic analysis, such as Wiener's Lemma, in order to 
study twisted convolution operators. Recently, in the authors described an new approach to classify 
the invertibility of £^-sequences with respect to the twisted convolution for rational parameters. 

In this manuscript we extend the idea of in the sense that we take a different approach which 
allows far better insights into the problem. Specifically, we only deal with sequences and show explicitly 
how efficient inversion schemes can be derived by rather simple (though sophisticated) manipulations of the 
twisted convolution. The essential idea is to split up the twisted convolution into a finite number of sums 
that can be incorporated into a special matrix algebra. In this matrix algebra we then prove a special type 
of Wiener Lemma which is the most challenging part from a mathematical perspective. 

The paper is organized as follows. The first section briefly outlines the basic definition of the twisted 
convolution. In this section we further discuss the example of twisted convolution on the finite group Zp x Zp. 
This example serves the purpose to motivate the introduction of the matrix algebra that appears in Section 

3 where we prove Wiener's Lemma for a special subalgebra. Section 4 links the twisted convolution to 
time-frequency operators. More specifically, it shows how the results shown in Section 3 can be used in the 
context of Gabor frames. In the last section, we give a short outline of the application of the presented 
approach for inverting frame-like Gabor operators. 

2 Twisted Convolution 

Let p and q be integers and relatively prime. We define the twisted convolution for sequences a,b £^{I?'^) 
by 

{a\\b)m,n= Yl ak,lbm-k,n-l^^'^~^^'^ (1) 

where uj = e^'^'^il'P and • denotes the inner product in M'^. Although the twisted convolution depends on 
p, q we do not specify this dependence because p, q will always be given and fixed beforehand. In Section 

4 we show how the twisted convolution is related to a class of operators with a special time-frequency 
representation. 
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In contrast to the conventional convolution with symbol *, in which to = 1, the twisted convolution is 
not commutative, and turns £^{7?'^) into a non-commutative algebra with the delta-sequence 6 as its unit 
element. 

We tackle the problem to study the invertibility of twisted convolution operators. Non-commutativity 
is the main subtle point in this problem. In fact, the question when the mapping 



for some 6 E ^ is invertible and how we can compute the inverse is more difficult than for a commutative 
setting. In particular, Wiener's Lemma which deals with the problem that if, for some b € £^ , Cb is invertible 
on £^ then the inverse is generated from an element again in i^, has to be proven separately. An abstract 
and more general proof of Wiener's Lemma for twisted convolution is given in ^Uj. Herein, we focus on a 
constructive method for studying the invertibility of the twisted convolution with the rational parameter 

q/p- 

In the following subsections we study the twisted convolution in a finite setting and draw analogies for 
approaching the problem of invertibility of Cb in the general case. 

2.1 Twisted convolution on Zp x Zp 

In what follows we describe the twisted convolution on the finite group F = Zp x Zp. The standard 
(commutative) convolution of two elements f,gG C^^p is defined by 



where operations on indices is performed modulo p. 

In analogy to the infinite case, we define the twisted convolution f \\g of two elements f,g(z C^^p by 



Cb-.aei 



p-i 

if * 9)m,n — ^ ^ fk,l9m—k,n—l 
k,l=0 



p-1 



(/tl5)m,n= ^ fk,l9m-k,n-l'^^"^ 



k,l=0 



with LO 



^2mq/p ^ For a fixed g, the twisted convolution can be seen as a linear mapping Cg : f ^ f \\g 



whose matrix G is block circulant with p blocks, i.e. 



Go Gp-i 



G — C{Go, Gp-i, . . . , Gi) 



Gi 



Go 



G2 



y Gp-2 ■ ■ ■ Go J 
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Each block has entries of the form 

iGj)ki = uj^''gj,k-i ■ 

Note that for the regular convolution each block is itself circulant. For the invertibility of block circulant 
matrices we apply a well known result from Fourier analysis. 

Lemma 2.1. 3 The matrix G = C{Gq, Gp-i, . . . , Gi) is invertible if and only if every 
s = 0, . . . ,p — 1, is invertible. In this case 

G ^ = G{Ho,Hp^i, . . . ,Hi) 

where Hr = ^ ^1=1 e'^^''''/P{Gs)-^ . 

By analyzing Gg, we see that all blocks are unitary equivalent, in the sense that 

where denotes the unitary matrix with entries 

1 if p — r = I — k, 
else . 

Since p and q are relatively prime, we obtain all blocks by such a unitary transformation. This implies that 
showing that if Go is invertible, then all Gg are invertible for s = 1, . . . ,p — 1. In other words, the p x p 
matrix Gq contains all the information about the invertibility of Cg. An easy computation shows that the 
entries of Gq are given by 

p-i 

(G'o)„,z = 5]a;"W-/- (2) 

A;=0 

We will later see that this observation motivates the matrix algebra that we introduce to study the invert- 
ibility of the twisted convolution. 

Now, also all Gj^ satisfy the same unitary equivalence. It follows that we can read from Gq^ the element 

which inverts the twisted convolution f ^ f \\g, i.e., ^9 = 9^ 9^^ = 

The twisted convolution on Zp x Zp serves as analogy for modelling the twisted convolution for the 
continuous and the finite dimensional case. 

3 Main results 

Our aim is to find a way to describe those sequences that have an inverse in t]). To this end we 

divide the twisted convolution into a finite sum of weighted normal convolutions of sequences that have 
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disjoint support. We define such a sequence a'"''* by 

ak,i if {k,l) =p (r, s) . 



k,l 



(3) 
else, 



where r, s G Zp. Obviously, a^'** is supported on the coset (r +pZ'^) x (s and a = Ylr seZ'^ ^^'"^ ■ 

a sequence a having a coset support only for one index, e.g., on Z'^ X (s + pZ'^), we simply write a''*. We 
write =p for denoting the equivalence of integers modulo p. The idea of slitting a sequence into a sum of 
sequences supported on cosets has first been introduced by K. Grochenig and W. Kozek in 

Lemma 3.1. Leta,b,c be ml^{I?'^). 

(a) For r, s,u,v € Z^, a^'* * b^'^ is a sequence supported on the coset (n + r + pZ'^) x (u + s + pZ'^) . 

(b) If c = c '^ is invertible in {l^ , *), then c^^ is also supported on x pZ'^. 
Proof. Let a'''" , b''''" be sequences in i^iZ^"^) and k,l € Z^. Then 

= ^ ^ io,'^'^)t,wib'^'^)k+pZd-t,l+pZ'i-w 

m,n£Z^ {t,w)=p(m,n) 

Since a'''* is nonzero only for {t, w) =p (r, s), and 6"'^ for {k — t,l — w) =p {u, v), we obtain that {k, I) has to 
be equivalent to {u + r,v + s) modulo p for a'''* * 6"'^ to be nonzero. 
To show (b), let c = c'^ be invertible and e be its inverse. Then 



S = c*e = c'° * ( ^ e''"*) = ^ 



where, by previous calculations, c '*' * e '** is a sequence supported on Z'^ x (s -\- pZ'^) for each s € Z^. Since 
6 = J2seZ'' ^ '^^ ^^'^ elements of the sum have disjoint supports, c'^ * e '* = 5''*. But since 6 '^ = for s 7^ 
and 6 '^ = 6, we conclude that 

6 s = 

c'°*e'' = < 

s / 

and therefore e = e''^. □ 
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With Definition ©, we obtain for u,v & Zp 

= ^ ^ 0,k,lK+pZd-k,v+pZd-l^^'^ * 

/ , (a )kA" )u+pi'i-k,v+pi'i-i^ 



T,s&Z^ k,i 

Ef , ru—r,v—s\ , r) • 

r,sGZ^ 



In a more compact notation we have 



(a 1^6)"'''= ^ a'^'**5"-^'^-"a;("-^)-". (4) 

We observe now that the upper indices in @ behave hke a twisted convolution in X Z^. What changes 
is that we have sequences as elements and standard convolution instead of multiplication. 

Motivated by the block circulant structure of the twisted convolution on Zp x Zp as described in the 
previous section, we introduce a new matrix algebra which is isomorphic to {£'^ [Z,'^'^) , \\). 

Before we do so, we fix an ordering of the elements from Zp. Let N = p'^ and X = {1, . . . , N}. Then, to 
each i G 2 we assign an element ki from Zp and set ki = (0, . . . , 0). We will often write instead of ki. 

Let {M,®) be an algebra of p'^ x p'^-matrices whose entries are ^^-sequences and multiplication of two 
elements A, B & A4 is given by 



The identity element Id is a matrix with 5 sequences on the diagonal. 



Theorem 3.2. Let 



Mo = [AeM Aij = ^ uj""- kja^'M-k, ^^^f andij 



Then Mq is a subalgebra of A4. 

Proof. Define a mapping (j): {£^,\]) —>■ {A4, ®) by 



{a))ij = uj'^-^^a'^^^^-^K (5) 
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Then 4> is linear, {(j){6))ij = J2m€Z'i'^"^'^^^^''^* = 5 li i = j and zero otherwise. So (j){S) = Id. We 
emphasize that the mapping (p has been motivated by the matrix Gq described in the previous section. For 

m,l,s€Z^ 
m,l,seZd 
m,l,seZd 

= ^ {^J-'^^a'-^^'''^) * ( ^ ^s-{m-l)^m-l,ki-s^ 

s&Z'i leZd m&Zd 



= ^( ^ ^k„ ■ m^m,ki-k„^ ^ ■ kj gl,k„-kj^ 

= 5^<^(%n * HO')n,j = {4>{b) ® Ha)). .. 



Therefore (p is an anti-homomorphism, that is, 

(t>{a\\b) = (j){b) ® (p{a). 

Hence Aio is an algebra, being an image of an anti-homomorphism. □ 

Before stating the main theorem, we explore properties of elements of A^o- For i,j G X and a matrix 
A G A^o we define a new matrix A{j, i) obtained from A by substituting the jth row of A with a vector of 
zeros having S on the ith position, and the ith column with a column of zeros having 6 on the jth position. 

Lemma 3.3. Let A G Mo- Then 

(a) det(74) is a sequence supported on Z'^ x pZ*^. 

(b) det(A(l,i)) is a sequence supported on 7,'^ x (fej +p7,'^) for i ^X. 
Proof. Let Sn be the group of permutations of the set T. Then 

N N 

det(A) = ^(-l)-J]^,(,)^,= ^(-irn( E uj^'-'^'a^^M^-^^) 
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Since a is a permutation of T, 

{ka{l) - kl) + {ka{2) - ^2) H h {k„{N) " ^Af) = 0. 

Therefore, by Lemma [3.11 Gmi,...,mjv is a sequence supported on the coset (X^iei™"* + pZ*^) x pZ*^. Since 
^jgjmj runs over all Z^, we see that det(A) is supported on the coset Z'^ x pZ*^, i.e., det(A) = det(^)'''^. 

In order to compute the support of det(^(l,i)) for i €l, let Sn-i denote the group of permutations of 
{2, . . . , N}. Then for i = 1, . . . , iV, 

det(yl(l,z)) = (-1)'^+^ (-l)''^<x(2),i * • • • * ^<xw,i-i * ^<x{m),*+i * • • • * ^<x{w),iv 

= ( — 1)'^ ^;"T-2-feiH hm.i ■ A;i_i+mi+i ■ fci+iH hmjv ■ fcjv ^ 



Since o" is a permutation of {2, . . . , N}, 

{ka{2) -kl)-\ h {k„(i) - h^i) + {k„(i+i) - h+i) H h (fc^(Ar) - /cat) 

= (^cr(2) H ^ ^<7(Ar)) - (fci + A:2 H h %) + fci 

= (^<7{2) H ^ ^(7(Af)) - (^2 H ^kN) + ki = ki. 

Therefore, by Lemma ITTl Gm2,...,mjv is supported on (^^2 +p'L'^) x {ki +pIJ^), and since each nii runs 
over all Z^, det(A(l,i)) is supported on Z'^ x {ki+pZ'^). That is, det(A(l,i)) = det(A(l, □ 

Now we are in the position to state and prove the main result 

Theorem 3.4. [Wiener's Lemma for M^q] Let A G M.q. If A is invertible in A4, then B = A^^ € M.q. 

Proof. Since A € is invertible, det(A) is an invertible sequence in and there exists a matrix 

B £ Ai such that A ® B = Id. By Lemma 13. 3| det(j4) = det(A) ''' and by Lemma 13.11 its inverse, 
e = det(A)^^, is also supported on the same coset, hence e = e''^. By Cramer's rule the inverse of A is given 
by 

Bij = det(^(j,i)) * e. 

We see that by Lemma 13.31 (b). Bi^i is a sequence supported on Z"^ X {ki +pZ'^). Let 6 be a sequence defined 
by 

b = Bi i + i?2,i + . . . + Bj\f^i. 



Then Sj^i = '^j^x^^^'^^ ■ Define a new matrix, denoted by B, as 

Then B G A^o and we will show that B = B, that is, B is the inverse of A. 
Since B is the inverse of A, 

Idi,i = {A®B).^=Y,Ai,j*hi 

j€l meZd; 



= ^ ^ ^ ^m,ki-kj ^ ^n,kj 

_ ^ ^ ^ ^(m-n)-kj^m-n,ki-kj ^ ^n,kj 

mez^ jeT nez^ 
= ^ G{m,ki), 

meZd 

where G{m,ki) = ^ pj ^„p^d w'-'"""'-' ' *^-''a"*~"'''^'~'^^ *?>"-'*^j is a sequence supported on (m + pZ) x (/cj+pZ). 
Therefore, G{ki,ki) = 5 and G(m,ki) = for m 7^ fci and i / 1. Using the above identity we will show 
that A ® B = ld, and by the uniqueness of the inverse we will conclude that B = B: 



{A®B).. = Y,As*Bs,j = 

s€I meZd; neZd; 

Tfl ■ ks f^^Tl ' kj ^TTl^ki ks ^ l^i^s kj 



meZ^ s€I neZ^ 



= ^ ^ ^ i^{ks+kj)-mjn-kj^m,ki-kj-ks ^yn,ks 



m&Z^ sex ngZd 



_ ^(ks+kj) ■ {m-n) ■ kj ^m-n,{ki-kj)-ks ^ yfi,ks 



mez^ sei nezd; 



— Lij'"''^'^^^ ^(m-n)-ks^m-n,{ki-kj)-ks ^yn,k 



6 ki — kj = ki <^ i = j; 

ki - kj ^ ki ^ i ^ j; 

Hence, A ® B = I. □ 



^ u}"'-^JG{m,ki-kj) = < 



Theorem 13.41 provides the key result to study invertibihty of twisted convolution. Indeed, for a given 
sequence a in £^ we look at the corresponding matrix A = 4>{a) as defined in ((21). If A is invertible in (A^, ®), 
which can be checked showing that the determinant is invertible in (^^,*), then its inverse A^"^ is of the 
form for another element 6 in This element 6, in turn, provides the inverse of a in {l^, t]). 

The approach is constructive in the sense that algebraic methods such as Cramer's Rule can be applied 
to find the inverse of A. Then, the sequence b can simply be read from the entries of A~^ according to 
the mapping (p. In particular for small p and d this method leads to fast inversion schemes for the twisted 
convolution operator. In the last section we will show explicitly how this works in the case of d = 1. 

4 Twisted convolution and Gabor analysis 

Central objects in time frequency analysis are modulation and translation operators. Although most of the 
upcoming notation can be given in the more general setting of locally compact Abelian groups we restrict 
ourselves to M*^ in order to simplify the readability of this article. 

For £ M'^ we define the translation operator and the modulation operator on L^(R°') by 

= /(--x), 
M^/(.) = e2™7(.), 

respectively. Many technical details in time-frequency analysis are linked to the commutation law of the 
translation and modulation operator, namely, 

M^T, = e2™-r,M^. (6) 

The time- frequency shift for j;, u; € M"^ is denoted by 

7r(x, w) = T^M^. 

It follows from © that 

This shows that time-frequency shifts almost allow a group structure. Incorporating the additional phase 
factor into a more extended group law leads to the so-called Heissenberg group. For more details about this 
topic, the reader is referred to [7]. 
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Gabor analysis deals with the problem of decomposing and reconstructing signals according to a special 
basis system which consists of regular time- frequency shifts of a single so-called window function jSj • Let 
A be a time-frequency lattice, i.e., a discrete subgroup of the time-frequency plane M?'^, and let g be in 
L'^{R'^). Then we define a Gabor system Q{g,A) by 

g{g,A) = {^(A)5|ag a}. 

We associate with this Gabor system the positive operator 

5:/GL2^5/ = ^(/,7r(A)5)vr(A)<7. 

AeA 

'2/ 



If the operator S is bounded and invertible on L (M ), then Q(g,A) is called a frame and S the associated 
frame operator, cf. [V . 

Many studies in Gabor analysis are devoted to the frame operator 8_. In what follows we will describe 
the so-called Janssen representation of such operators. To this end we need the notion of the adjoint lattice, 
i.e., 

A° = { A° € M^^ I ^(A)^(A°) = ^(A°)7r(A), A G A } . 
In Pl 1^ ^2 it is shown that the frame operator S satisfies Janssen representation, 

S= (5,MA°)5MA°). (8) 

A°eA° 

At this point, the question arises if we can deduce the invertibility of the operator S from the Janssen 
coefficients {{g, iT{X°)g)). It is known from frame theory that if S is invertible, then its inverse is of the same 
type, that is, it also has a Janssen representation. 

In order to better understand the main ingredients of this problem we transfer the model to an operator 
algebra. To this end we restrict our discussion to so-called separable lattices of the form 

A = aZ'^ X 

for some fixed positive numbers a and (3. An easy computation based on ((T)) shows that 

A° = p-^IJ^ X a-^U^ . 
We define the operator algebra A as in by 

A= [S= akMP~'k,a-H) I a = (a,.,;) G l^Z^'') } . 
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The restriction to £ -sequences guarantees absolute convergence of the sum of time-frequency shifts. Let k 
be the mapping 

K : a £ ^ K{a) = ak^iTT{f)^^k,a^^l) £ A. 

k,l&d. 

Then, as aheady observed in ^T], we have 

n(a)K(b) = n{a\\b) 

and k{6) = Id where S and Id denote the Dirac sequence and the identity operator, respectively. Both 
represent the unit element of the corresponding algebra. It follows that k is an algebra homomorphism from 
t]) to A, and invertibility of an element in A can be transferred to the invertibility of the associated 
^"'^-sequence with respect to the twisted convolution. 

It is important to observe, that all the results go through also for weighted £^-spaces. These facts are 
used to design dual Gabor windows of a special type, cf. [TU] . 

In the following section we give an example of how this approach can be explicitly used in Gabor analysis 
of one-dimensional signals. 

Remark. The above results, with the help of metaplectic operators, carry over to the more general class of 
lattices, called symplectic lattices. A lattice C M^"' is called symplectic, if one can write A^ = DA where 
A is a separable lattice and T) € GL2d{^)- To every T> € GL2(i(M), there corresponds a unitary operator 
/x(!D), called metaplectic, acting on L'^{W^). One can show that a Gabor system on a symplectic lattice is 
unitary equivalent to a Gabor system on a separable lattice under niV), and 

Hence, to analyze the invertibility of a frame operator S associated to the window function g € L^(M'^) and 
symplectic lattice A^, it suffices to analyze a frame operator associated to ^{T>)g and a separable lattice A. 
For more details see |H]. 

5 Application to one-dimensional signal space 

In this section, we briefly describe how the presented inversion scheme applies to Gabor frame operators in 
a one-dimensional setting. A more detailed discussion also for flnite dimensional signals is described in ^H] . 
Assume d = 1. Let a be in and a, j3 be constants such that aj3 = p/q with p, q relative prime. Set 

k,iez 
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In order to verify if K{a) is invertible on L^(R) we simply look at the coefficient sequence a and check 
whether a is invertible in {£^(7?), t]). To this end, we apply the above results and switch to the matrix A 
whose entries are defined by 

p-i 

m=0 

with Lo = e^'^*^/^. Next, we need to show that the matrix A is invertible in (A^,®). For example, we can 
calculate the determinate which is a sequence in and show that it is invertible in (i^,*). 

Assume that the determinant of A is invertible. We denote its inverse by e. By Cramer's Rule, we 
compute the first column of the inverse matrix S of A as 

Bk,o = det ^(0, k) * e, 

ioT k = 0,...,p-l. Then 

p-i 

k=0 

provides the inverse sequence of a which, in turns, gives K{a)~^ = 

Note that for p = 1, the twisted convolution turn into normal convolution and we can simply apply the 
standard Fourier inversion scheme of sequences in (£^(Z^), *) since in this case the matrix A reduces to the 
sequence a. 
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